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Abstract 

I present a systematic study of self- assembled InAs/InP and InAs/GaAs quantum dots single 
particle and many body properties as a function of quantum dot-surrounding matrix valence band 
offset. I use an atomistic, empirical tight-binding approach and perform numerically demanding 
calculations for half-million atom nanosystems. I demonstrate that the overall confinement in 
quantum dots is a nontrivial interplay of two key factors: strain effects and the valence band offset. 
I show that strain effects determine both the peculiar structure of confined hole states of lens type 
InAs/GaAs quantum dots and the characteristic ,, shell-like" structure of confined holes states in 
commonly considered "low-strain" lens type InAs/InP quantum dot. I also demonstrate that strain 
leads to single band-hke behavior of hole states of disc type (,, indium flushed") InAs/GaAs and 
InAs/InP quantum dots. I show how strain and valence band offset affect quantum dot many- 
body properties: the excitonic fine structure, an important factor for efficient entangled photon 
pair generation, and the biexciton and charged excitons binding energies. 
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I. INTRODUCTION 



Fully ab-initio, parameters free, modeling of million atom self-assembled^ or nanowire^ 
quantum dots (QDs) is still beyond the reach of current computers. For practical, atom- 
istic calculation, semi-empirical approaches like the empirical tight-binding (ETB)^"— or 
the empirical pseudopotential method (EPM)^~— are typically employed. The calculation 
scheme usually starts with strain field calculation followed by the single particle calculation 
followed then with the configuration interaction approach to obtain many-body (exciton, 
charged exciton, multi-exciton) spectra^^^^^. 

Semi-empirical approaches use sets of fitted parameters determined to reproduce bulk 
properties like effective masses, bulk deformation potentials and gaps at diflFerent points of 
the Brillouin zone^^. Bulk derived parameters are later used for calculation of nanosize sys- 
tems. Apart from potential (,,bulk to nanosystem") transferability issues, one may question 
the reliability of important parameters describing bulk electronic structure, that act as the 
input data for the empirical fitting procedure. For example, in the case of the absolute 
In As valence band deformation potential (a^) not even the sign of this quantity is unam- 
biguously determined—"—. As semiconductor (self-assembled or nanowire) quantum dots 
are typically mixed material systems an additional empirical bulk parameter, the valence 
band oflFset between quantum dot and surrounding matrix material, has to be incorporated 
into the Hamiltonian. This ,, natural"— ^ valence band offset (VBO) determines the depth 
of unstrained hole and electron confining potentials and combined with strain and defor- 
mation potentials constitutes the overall confining potential for the strained case, i.e. the 
,, strained" band offset. There is again a substantial uncertainty of the natural VBO values, 
e.g., the reported InAs/GaAs VBO varies from 50 meV to 500 meV— The source of 
this discrepancy is not only due to the difference between experimentally and theoretically 
reported values, but^: "It is to be emphasized, however, that even within the framework 
of Kohn-Sham DFT, different computational schemes result in different predictions for the 
natural band offsets [...]". Y.H. Li and co-workers^^ state additionally that there has been 
" long-standing anomalies between theory and experiment" and in their ab-initio calculation 
"For GaAs/InAs the predicted offset is increased from 0.06 eV in the previous calculation 
to 0.50 eV". 

For a given pair of materials (e.g. InAs/GaAs) the VBO has only one value, however 
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as discussed above, this value is practically unknown or given with large (> 100 meV) 
uncertainty. E.g. the VBO=0.06 eV is used in the EPM approach, the VBO?^ 0.17 eV is 
utilized in k-p studies^^, the VBO=0.23 eV is incorporated in the ETB model by Boykin et 
al."^^, while VBO?^ 500 meV is reported by the DFT calculations^^. In this work, as practical 
resolution of this problem, I utilize an approach in which I perform calculations by an 
artificial variation of the VBO over wide range of values. While the VBO has been considered 
as a merely technological parameter of lesser importance, the current study shows that it 
is quite the opposite. Both the single particle and the many body properties are aflFected 
by the choice of the VBO and the caution should be exercised before applying different 
VBO values in the semi empirical calculation. Abstracting form the experimental reality 
(there is only one, yet unknown VBO value, for a given pair of materials), the "artificial" 
modification of the VBO is by itself a very interesting theoretical tool to study the effects of 
band confinement versus other effects such as strain. For example, the calculation of spectra 
for InAs/InP and InAs/GaAs lens type quantum dots using the same (yet artificial) VBO 
helps to understand the difference between both types of nanosystems. 

The ambiguity of the ay and VBO bulk material values may thus affeclp^ the accuracy of 
the qualitative description of confined valence band states in semiconductor nanosystems. 
This is further important as recent empirical pseudopotential method works (e.g. Refs.— 
or— etc.) still refer to older EPM parameterizations^^ utilizing the " questioned" VBO 
valuer 50 meV. The Section E of the paper discusses this point showing that while 'natural' 
VBO in ETB and EPM can differ significantly, the strained band offsets in two approaches 
are very similar. 

Despite over a decade of intensive studies the problem of detailed understanding of holes 
in self- assembled InAs quantum dots is still an active field of research^^. Ediger et al.— 
observe characteristic spectral structure of hole states in InAs/GaAs lens type quantum 
dots leading to a non trivial hole charging pattern of excitonic complexes. In our earlier 
work^^ we have noticed that InAs/GaAs quantum dot shape affects spectral properties of 
holes significantly. Recently Gong et al.^^'^^ have calculated electronic structure of InAs/InP 
lens type quantum dot and speculated that differences with respect to analogous InAs/GaAs 
system is due to different common ion or different VBO in both systems. In this work, by 
the systematic VBO analysis, I show it is not the VBO, but rather strain effects that play a 
dominant role in determining the character of single particle holes states in quantum dots. 
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A two photon cascade from the QD biexciton state can generate entangled photons^'' and 
has attracted great interest for apphcations in quantum information. However, anisotropic 
exchange sphtting of quantum dot bright excitons, induced by the asymmetry of typical 
quantum dots, inhibits entanglement. Understanding the origins of fine structure splitting 
is thus of great importance for potential quantum dot applications. Recently, yet another 
scheme for entangled photons generation has been proposed^^i"— , in which the biexciton 
binding energy is tuned to zero. As in a typical experiment different excitonic complexes 
(both charged and neutral) are observed together, the prediction of the spectral line order— 
or binding energies is usually far from trivial. To solve such issues, an inverse approach for 
quantum dot calculation has been recently proposed^. In this method one uses excitonic 
spectroscopy experimental data to determine quantum dot structural properties. However 
the accuracy of such prediction must depend on the accuracy of the many-body calculations 
and indirectly on empirical parameters (such a the VBO) used in the calculation. 

In this work I compare properties of strained and unstrained systems, study the role 
of quantum dot shape and the evolution of single particle and charge probability densities 
as a function of valence band offset. By a systematic VBO analysis, I show that strain 
and valence band offset effects play different, important roles in determining the character 
of single particle states in quantum dots. Finally, I show that the choice of VBO affects 
substantially many body energies, in particular biexciton and trions binding energies and 
the excitonic fine structure. 

II. SYSTEMS AND METHODS 

In the following I present a systematic study of lens and disc type In As quantum dots, 
surrounded by either InP or GaAs, as a function of valence band offset, with strain effects 
included or artificially neglected. The height of the disc type quantum dot is /i = 3 nm and 
the base diameter is = 16.8 nm. For the sake of comparison with the EPM calculations2i"i2i 
the height of the lens type dot is chosen as /i = 3.5 nm and the base diameter is = 25 nm. 
Both dots are located on 1 nm thick wetting layer. The presence of the wetting layer 
is particularly important for disc type quantum dots as it lowers the overall quantum dot 
symmetry from to (the lack of " rotoinversion operation" and therefore both kinds 
(lens and disc) of quantum dots have low symmetry. 
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The calculation consists of several major steps: first atomic positions (strained and un- 
strained) are calculated. There is a lattice mismatch between dot material (InAs) and 
surrounding matrix material (GaAs or InP). To calculate strain relaxed positions I use the 
atomistic valence force field (VFF) approach of Keating^^. This method is described in more 
detail in ref.— and in our previous worksl^^i^. The size of the computational domain, in- 
cluding more than 50 millon atoms, guarantees convergence of the strain distribution^-. 

Due to small lattice mismatch of InAs and InP I neglect the piezoelectric effects in the 
present calculation, following similar arguments by Gong et al.^^ who ignore piezoelectricity 
in the empirical pseudopotential work on InAs/InP quantum dots. Consistently, piezoelec- 
tric effects can also be neglected for low aspect ratio^*^*^ lens and disc type InAs/GaAs 
quantum dots, where the piezoelectricity is either negligible^^ or the contribution due to 
second-order effects tends to cancel linear terms^^^. In particular for lens type InAs/GaAs 
ref.— states that: "for smaller aspect ratios [...] first- and second-order effects compensate 
each other with respect to their impact on the electronic states." Finally the piezoelectric- 
ity is neglected for a sake of a fair comparison with other approaches where this effect is 
neglected (comment 41 from ref.^^). 

In the second step of the calculation, the single particle states are obtained by building 
the sp^d^s* tight-binding Hamiltoniani^^^ and then diagonalizing the Hamiltonian by means 
of the Arnoldi algorithm with the matrix- vector multiplication parallelized using OpenMP 
approach on 48 core, shared memory system. For the (artificially) unstrained^^ systems I 
use bulk TB parameters and thus there is no strain contribution in the Hamiltonian nor the 
relaxation of atomic position is accounted for. For the strained systems, since strain effects 
change bond lengths and angles, strain relaxed positions are used to modify TB parameters 
(diagonal and off-diagonal matrix elements) following the description given in detail in my 
earlier work—. 

I have calculated the single particle spectra for many VBO values from 50 meV to 500 
meV with 10 meV step. In this work I use a multi-scale (multi-domain) approach where 
smaller computational domain is used for the single particle calculation^*^. Yet, as the 
number of atoms in the TB domain is larger than half a million, in order to make the 
entire process feasible I calculate only eigenenergies of several lowest electron and holes 
confined states. Then for several chosen VBO values I additionally calculate eigenstates and 
plot corresponding probability charge densities. Finally, for several VBO values electron 
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and hole Coulomb matrix elements (Coulomb direct and exchange integrals) are calculated 
according to the approach given in ref.^^ and the many-body Hamiltonian for the exciton 
(X), the biexciton (XX), positively (X~)and negatively charged (X^) trions is solved using 
the configuration interaction approach^'^. 

The quantum dot and the surrounding matrix may share the same anion (e.g. InAs/GaAs 
QD), cation (e.g. InAs/InP QD) or not have a common ion (e.g. InAs/GaP). In the empirical 
tight-binding the treatment of quantum dot/material interface atoms is ambiguous. In 
Boykin et al.^^ approach, this is handled during the fitting procedure where the diagonal 
Hamiltonian matrix elements of the common atom are kept the same in both materials. 
The value on the material band offset is incorporated into off-diagonal hopping matrix 
elements. This approach removes the necessity of modifying on-site matrix elements for 
interface atoms, but complicates significantly for ternary systems like InAs/InP/GaP, where, 
e.g., the bulk GaP properties are indirectly coupled to InAs bulk properties through the 
fitting process. Additionally as VBOs are embedded into the tight-binding parameterization, 
a necessary, complicated refit would be needed for every different VBO values. In this work 
I use the approach^^ in which I account for the valence band offset by shifting diagonal 
matrix elements of the quantum dot material. On-site matrix elements of interface atoms 
are calculated as a weighed sums of neighboring atoms materials on-site parameters. The 
Hamiltonian is built and then diagonalized for each VBO value: this gives a capability of 
freely tuning VBO, suitable thus for the VBO dependence studies. 

III. RESULTS 

Figure [1] shows the evolution of single particle electron and hole energy levels for the lens 
type InAs/GaAs quantum dot as a function of the InAs/GaAs VBO. With the increasing 
VBO value there is an energy upshift due to overall ,, reference level" (VBO) energy shift, 
while the effective band gap Egap = ei — hi does not change considerablj^^*^. In this work 
I show that while the above statement is generally true^^i^*^, the situation is far more 
complicated for spectral quantities other Egap. To analyze the specifies of electron and hole 
spectra in the following plots I subtract the corresponding carrier ground state energy as 
illustrated on the example Figure [lb) and [lb) . 

Valence band offset is related to (natural or unstrained) conduction band offset CBO 
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through the foUowing relation CBO = matrix gap — dot gap — VBO^ where matrix gap is 
the bulk band gap of the matrix material (GaAs or InP) and dotgap is bulk band gap of 
the quantum dot (InAs) material. For the VBO in 50 — 500meV range, the CBO varies 
correspondingly from ^ leV to ^ 0.5eV, with little difference between InAs/GaAs and 
InAs/InP cases, due to 0.1 eV difference of InP and GaAs bulk band gaps. Increasing the 
VBO corresponds thus to decreasing of the CBO, i.e. lower confinement of electron states. 
Even though CBO is nominally larger than VBO, effects of CBO variation on electron states 
should be observable due to their lower confinement (effective mass). 

Figure [2] and [3] show energy levels corresponding to several lowest electron (upper-(75) 
and hole (lower- 1/5) states calculated for InAs/InP and InAs/GaAs disc type and lens 
type quantum dots as function of quantum dot-matrix valence-band offset (VBO) with 
strain effects either artificially neglected or included. Figures from Figure d] to Figure [TOl 
show corresponding charge/probability density isosurfaces. These figures contain substantial 
amount of information and will be analyzed in detail in the following part of the text. 

A. Electron states - strain effects neglected 

In an artificially unstrained InAs/(InP,GaAs) disc type quantum dot energy spectra of 
lowest electron levels reveal shell like structure [Fig. [2] a) and b)], with the ground electron 
state of s character [Fig. [4]. Despite the absence of strain and the rotational shape symmetry 
of the disc quantum dot, the presence of atomistic interfaces and low symmetry of underlying 
crystal lattice introduces the asymmetry into the system Hamiltonian^^. Thus, there are 
two closely spaced (splitting < 1.5meV) excited states (e2 and 63) of p-like character (of 
approximate angular momentum character L = ±1), followed by two excited (64 and 65) 
closely spaced states (splitting ^ 4 meV) of d-like character. Splittings within p and d shell 
are still however much smaller then spacings between different shells: s — p (60 meV) and 
p — d (70 meV). The (of "2^" character^^) is separated from the lower lying rf-like states 
64,65 by ^ 20 meV, a hallmark of disc-like confinement^. Such structure of levels with 
quasi-degenerate energies corresponding to L = 0, ±1, ±2, ... is to be expected for nominally 
cylindrical disc-shaped quantum dots. 

Charge probability densities corresponding to several lowest electron states in artificially 
unstrained InAs disc type quantum dot practically do not change in the considered range 
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of VBO values and are very similar for both InP and GaAs matrices [FigSj. For small 
VBO ^ 100 meV there is a slight elongation of the p-shell states (e2 and 63) along [110] and 
[110] crystal axis, but otherwise these states have well defined cylindrical-like symmetry. 
For the VBO changing from 0.1 eV to 0.4 eV, the corresponding ground electron state 
localization inside quantum dot drops only by 3%, i.e. from 81% to 78% in InAs/InP case 
and from 86% to 83% in InAs/GaAs case. With no strain included, the confinement is 
generally somewhat lower for InAs/InP systems due to lower InP band gap and thus lower 
CBO than for the InAs/GaAs case. 

For the disc type quantum dot, in the absence of strain, spacings between diflFerent shells 
s — p and p — d do not change much as a function of VBO, however Figure [5] shows that the 
splitting of the p-shell (,, p-shell anisotropy") increases monotonically (quasi-parabolically) 
with decreasing confinement, most likely due to the increasing role of interface effects in 
progressively shallower (decreasing CBO) well for electrons. Notably, in the absence of 
strain, the p-shell splitting is systematically larger in disc type InAs/InP quantum dots 
then in disc type InAs/GaAs systems, which I speculate, can also be related to 5%) 
lower confinement of electrons in InAs/InP quantum dots. 

In the effective mass approximation, lens type quantum dots are expected to show 2D 
harmonic oscillator-like spectrum^'^^. In my atomistic calculations for the lens type quantum 
dot [Fig. [3] a) and b)] I also observe the well pronounced shell structure of electrons levels 
and well defined nodal-structure of corresponding charge densities (Fig. [6]). For the lens type 
quantum dot, electron s — p and p — d level spacings slightly decrease with the increasing 
VBO. The ground electron state charge distribution is apparently not affected by the choice 
of the matrix material. The lower electronic p-shell (62) state is localized along [110] crystal 
axes and the higher p-shell state is localized along [110] axes for the InAs/GaAs system, while 
for the InAs/InP system both p-shell states maintain cylindrical-like symmetry. Figure [5] 
shows also that with strain effects neglected the p-shell splitting depends more strongly on 
the absolute depth of the confining potential due to the VBO rather than on the particular 
quantum dot shape. Higher lying states are approximately of the same symmetry for both 
GaAs and InP matrices (Fig. [6j) with similar inter-shell spacings 60 meV). 
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B. Electron states - strain effects included 

With the strained effects accounted for, the first important difference is a significant 
increase of the p-shell splitting for both disc and lens system (Figure [5j) and the well pro- 
nounced elongation of the p-shell states along [110] and [110] crystals axes (Figured]). Inter- 
estingly the p-shell states elongation for the disc type quantum dots is opposite to that of the 
lens type quantum dots (Figure [6]). Also, for the lens type quantum dot the orientation of 
the p-shell states is reversed^S when compared to the strain-free case: the anisotropy due to 
atomic interface is thus reversed by the anisotropy due to strain. For the disc type quantum 
dot, the competition of two different anisotropy sources manifest itself by non-monotonic 
change of the p-shell splitting as a function of the VBO (Figure [5]), especially pronounced for 
the highly strained InAs/GaAs system. With strain effects included the p-shell splitting is 
generally higher for lens type then for disc type quantum dots, most likely due to larger sur- 
face (interface) /volume ratio in lens type systems. With strain effects neglected the p-shell 
splitting was generally larger for InAs/InP systems, however with strain effects included this 
trend reverses, the splitting is dominated by strain, and therefore is larger for InAs/GaAs 
then for InAs/InP systems. To summarize, the overall structure of p-shell electron states 
is thus determined by combined: matrix material, dot shape and VBO value and none of 
these factors can be neglected^. 

For confined quantum dot states (both electrons and holes) spacings between different 
shells (s-p, p-d, etc.) increase with the increasing confining potential depth (i.e. band 
offset, the GEO and the VBO for electron and holes correspondingly) as expected from 
the quantum confinement effect [Fig. [2] and Fig. [3]. However splitting of levels within a 
given shell decreases with the increasing confinement. This is due the progressively larger 
localization and effectively smaller influence of the material interface, which acts as the 
source of splitting. Gonsequently, higher lying (d-shell) states properties are even more 
susceptible to the choice of the VBO due to their lower confinement. 

C. Hole states - strain effects neglected 

Interestingly even when neglecting the strain effects, the ground hole state of the consid- 
ered disc and lens type quantum dots has well defined s-like symmetry [Fig. [7] and Fig. [8] 
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and is predominately of heavy-hole character. This initially may come as a surprise as in 
such case there are no heavy hole-light hole splitting terms due strain in the Hamiltonian^. 
Yet, the quasi- two-dimensional confinement in fiat quantum dot systems is efiicient enough 
to separate both type of holes. Alternatively, one can associate heavy hole states with the 
in-plane component dominated by Px and Py atomic orbitals^^^^^ and affected by small, lat- 
eral confinement. Then the light-hole states are the predominantly constituted by orbitals 
and highly infiuenced by the vertical confinement, and thus energetically shifted away from 
the ground hole state. 

Higher lying states however, for both types of quantum dots, reveal complex, mixed 
angular momentum character [Fig. [3 and Fig. [8] and show no clear shell-like structure of 
their energy spectra [Fig. [Hand FiglGj, that was so characteristic for the single-band-like 
confinement of electron states. For lens type dots, the first and the second excited hole 
states are not even of the p-like symmetry and no nodal planes are observed. 

Due to strong (> 90%) confinement of hole states in quantum dot area, with strain 
effects neglected, there is little difference between hole states properties with respect to 
the surrounding matrix (GaAs/InP). This is particularly true for large VBO values, and 
better confined states, where the interface effects play a lesser role. And consequently there 
is almost one to one correspondence of charge densities for the disc type InAs/InP and 
InAs/GaAs quantum dot cases (Figure [71). 

Interestingly, even for the VBO= eV and no strain effects included I still observe 
quantum dot hole confined states due to material property (,, effective mass") discontinuity 

(Fig. ED. 

D. Hole states - strain effects included 

Strain affects holes states properties significantly. For disc type InAs/GaAs and InAs/InP 
quantum dots, in the " realistic" range of VBO values (210 — 350 meV). the structure of 
confined hole states is in a vivid contrast to the strain-free case and resembles that of single- 
band electron confined states, with well pronounced shells of p-like and d-like symmetry. 
This effect is well visible, both in the charge distributions [Fig. [9] and in the energy spectra 
[Fig. [2] g) and h)]. The hole p-shell splitting lies within few meV, i.e. much smaller then 
(25 — 30 meV) s-p level spacing. Strain lifts heavy hole-light hole degeneracy^^^^^^ and 
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thus effectively decouples light-hole component of the confined hole function, leading to the 
single band-like behavior of hole states of disc type (,, indium fiushed"^^) quantum dots and 
characteristic shell-structure know also from the experiments^. Similar conclusion may also 
be drawn for the disc-shaped nanowire quantum dots (which are however not placed on the 
wetting layer). While we have obtained this characteristic spectra in our earlier work—, 
we have attributed it the quantum dot shape. In this paper, I demonstrate that strain, 
rather then shape only, is responsible for the typical spectra of the disc type quantum dots. 
Additionally, it is actually the strain that also leads to the characteristic harmonic-oscillator- 
like structure of confined holes states in lens type InAs/InP quantum dot as seen on Fig. [3] 
g) and Fig. [101 where no shell-like structure could be observed with strain effects artificially 
neglected [Fig|3]e) and Fig|8]. Strain cannot be neglected even for InAs/InP quantum dots, 
commonly considered ,, low-strain" systems. 

Holes states properties may vary significantly in the considered range of VBO values, as 
VBO = means no confinement other then due to strain and material properties (,, effective 
mass") discontinuity (Fig. [9]). For small VBO values (and lesser confinement) there is a 
significant leakage of the hole wavefunction into the surrounding matrix and into the highly 
biaxially strained wetting layer [Fig. [9] and Fig. [10]. Strained disc type quantum dots reveal 
thus strong dependence on the choice of valence band offset, however for VBO > 300 meV 
their spectral properties stabilize. 

Higher lying holes states have mixed angular momentum character and their evolution 
with respect to the VBO seems to depend in a complicated way on the relative evolution 
of different angular momenta components. In terms of states localization I can label hole 
states by those localized predominately along one of two non-equivalent direction [110] and 
[110] axis correspondingly. These two species seem to evolve differently under the VBO 
change, leading to observed "levels crossings" [Fig. [2]g) and h)]. It is important to reiterate 
at this point that modification of the VBO is "artificial", yet it should be also noticed that 
is constitutes a very interesting theoretical tool. 

E. Strained valence band offset 

As expected, InAs/GaAs systems are affected by strain effects more than InAs/InP due to 
the large lattice mismatch of the former. For a small value of VBO = 50meV^ the confining 
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potential for holes is dominated by the strain contribution and the character (anisotropy) of 
this term leads to the ground (and the first excited) hole state of the apparent, yet unusual 
p-like symmetry (Fig. [TOj). I demonstrate results of the calculations for the particularly low 
(50 meV) VBO value as this number is customarily ,,used" in the empirical pseudopotential 
method (EPM^^^^^). However, I point that it is the strained band offset (being the ,, combina- 
tion" of both VBO and ay valence band deformation potential) that ,, enters" the calculation 
as the actual hole confining potential^^. While the ,, natural", i.e., strain free VBO and 
diflFer significantly between two methods (ETB: = 1 eV and EPM: = — 1 eV), impor- 
tantly, the strained band oflFset obtained by using two sets of and VBO parameters is 
very similar ^ 330meV^^. This is illustrated on Figure [TT] showing the confining potential 
profiles calculated for a InAs/GaAs lens type dot using the Bir-Pikus model"^^. This plot 
was obtained utilizing two sets of parameters: the ,, recommended" ones from the review 
paper by Vurgraftman et al.^^ (a^ = +1.0 eV, VBO = 0.21 eV) and those reported^^ by 
the EPM method (a^ = —1.0 eV, VBO = 0.05eV). Thus, in result, the effective confining 
potentials in the ETB and the EPM approaches are quite similar despite noticeably different 
bulk ("intermediate") target values. A more quantitative comparison of the quantum dot 
results obtained by the ETB and EPM was presented in our recent work-i^. 

The main difference between the lens type and disc type quantum dots comes from the 
fact that disc type quantum dots are subject to "smooth", slowly spatially varying strain 
(Fig. [T2l) . which also induces uniform heavy hole-light hole splitting leading to single particle 
like description of hole states. On the contrary, the lens type quantum dots are affected by 
spatially fiuctuating strain due to the curved shape of the quantum dots and the presence 
of jagged, step-like material interface (Fig [12]). Electron states are built predominately 
from s-type atomic orbitals and are affected only (in the Bir-Pikus formalism^^^i^) by the 
hydrostatic strain {Tr{e) = e^x + ^yy + ^zz)) that leads to simple energetic upward shift due 
to bond lengths deformation in the strained system. Thus electron states are not affected 
by highly spatially variable biaxial strain. 

On the contrary, heavy-holes confining potential (again in the Bir-Pikus formalism) 
along growth ([001]) quantum dot axis going through the dot center, is given as^-: 
Ehh = CLyTr (e) — bB (e), where B (e) = e^z — {^xx + ^yy) /2 is the biaxial component of 
strain, ay absolute valence band deformation potential and b is biaxial deformation poten- 
tial {bjnAs = —1.8 eV^^). Thus, for hole states, the strain related potential shift is an 
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interplay between the hydrostatic and biaxial strain, i.e. the biaxial strain {—bB (e) > 0) 
deepens the confinement for holes, the hydrostatic strain shallows the confining well in the 
dot region [a^Tr (e) < 0). 

FigJTSlshows the heavy-hole (in-plane) confining potential (obtained with bulk parameters 
from ref.^^) calculated for all quantum dots considered in this work. For the lens type 
InAs/GaAs quantum dot one can notice non-regular, "jittered" spatial dependence that will 
contribute to the splittings in the shell structure seen in the earlier part of the paper. These 
oscillations reach up to 50 meV and even lead to formation of attracting well for holes at the 
edge of quantum dot. This situation is further complicated at the interface as the biaxial 
strain changes sign^^ and reverses light-heavy hole ordering, thus affecting the valence band 
mixing. The interface related oscillation seems to play a small role for the relatively ,, low- 
strain" and high VBO value InAs/InP lens type quantum dots, yet it is the complicated 
character of strain that actually makes the lens type InAs/GaAs holes spectra so different 
from other quantum dots. 

This peculiarity leads, e.g., to the non-Aufbau hole charging pattern and been confirmed 
by the experiment-^'^^, the well-established EPM and also recent k.p calculations^^. It is only 
for quite large^^ InAs/GaAs VBO values (> 0.450 eV) that hole charge distributions of lens 
type quantum dot resemble [Fig. [10] that of the harmonic oscillator-like states. However even 
for VBO> 0.7 eV the hole p-states splitting does not drop below 5 meV. One can expect this 
effect should be more pronounced in tall lens, cones or pyramid type (high aspect ratic'')of 
quantum dots^. 

If accurate (order of meV's) modeling of hole states is important, then the choice of 
the ,, natural" VBO plays an important role. This is particularly noticeable for lens type 
InAs/GaAs quantum dot [Fig. [3j. In the "recommended" range of VBO values^^ (210 — 
350 meV), hole hi — /12 level spacing varies between 9 meV and 14 meV and spacings of 
higher lying level change even more substantially, e.g. /13 — /14 spacing varies from 1 meV to 
7 meV. 

IV. MANY-BODY STATES 

Strain effects and the valence band offset play a fundamental role for the single parti- 
cle states in quantum dots. Next, I calculate many-body properties of single exciton and 
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excitonic complexes. 



A. Exciton fine structure 

Fig. [m shows excitonic fine structure^^ calculated for different VBO values for lens and 
disc type InAs/GaAs and InAs/InP quantum dots. The energy difference between the two 
bright excitonic states, the so called bright exciton splitting (BES), is related to the confining 
potential anisotropy^^. The BES is thus larger for curved-shape, highly anisotropic, lens type 
quantum dots, reaching values varying between 40 — 60 meV, while for disc type quantum 
dots the BES does not exceed 20 meV. These results are in quantitative agreement with 
recent EPM calculations^^, however the EPM systematically predicts much lower BES value 
then those reported in the experiments^ or in current ETB work. 

For lens type quantum dots the BES varies noticeably with respect to the VBO value, 
while it has almost flat VBO dependence for disc shaped systems. While Fig. [14] is the 
result of the full many-body configuration interaction calculation^^, including single particle 
states up to the d-shell, the majority of the BES splitting and the overall trend is very 
well reproduced (not shown here) by a single electron-hole anisotropic exchange (complex) 



where (e) electron and (h) hole are in their ground s-states and arrows correspond car- 
rier quasi-spins. V^^^ is responsible for mixing two bright excitonic states (tJ| and ^-l^and 
therefore leads to the BES. The BES is thus a two-body effect of the s-shell electron and 
hole spatial anisotropy and as such is not directly correlated with the electron p-shell split- 
ting/anisotropy (Fig. [5]). 

The BES is larger for highly strained lens type InAs/GaAs quantum dot (Fig. [Hj), as the 
low symmetry {C2v) of the confining potential is clearly related to strain^^. On the other 
hand the electron-hole anisotropic exchange interaction is also related to the electron-hole 
"overlap" in V^^^ . Therefore, similar confinement of electron and hole may lead to larger 
BES. Thus, for low strain InAs/InP disc type quantum dot, where electron and hole charge 
distributions ("envelopes") are almost identical (as seen previously on FigHland Fig|9j), the 
BES value (?^ 18/ieV) is almost two times larger than for highly strain InAs/GaAs disc type 



integral- 
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quantum dot 8/xeV). This is in clear contradiction to the recent EPM calculations^^, which 
predict highly reduced fine-structure splitting in InAs/InP quantum dots. To summarize, 
the BES is a non-trivial function of quantum dot shape, shape and quantum dot/matrix 
materials. 

The dark-bright exchange splitting is determined predominantly by a (real) exchange 
matrix element which also conserves spin (V^^~^). Bright-dark exciton splitting does not 
varies significantly with VBO: increased confinement of hole seems to be assisted by the 
decreased confinement of electron, and the overall electron-hole exchange integral and dark- 
bright exciton splitting does not vary too much. Based on similar arguments as for the BES, 
bright-dark exciton splitting for disc type quantum dots is generally larger. Additionally 
bright-dark exciton splitting will be larger for those quantum dots with better electron-hole 
overlap: InAs/GaAs rather then InAs/InP disc type quantum dot and InAs/InP rather then 
InAs/GaAs lens type quantum dot. 

Finally, dark exciton splitting is a pronounced and complicated function of the VBO, 
especially for InAs/GaAs systems, varying over an order of magnitude for the studied VBOs. 
Dark exciton splitting is generally larger for smaller VBO value, however no clear trend 
regarding the matrix material is observed. 

B. XX, X~ and X+ binding energies 

Fig. [T5l shows biexciton (XX) and charged trions (X~ and X+) binding energies, measured 
with respect to the single exciton energy {E§^ = E^^ — ^ E§ = E^ — E^ ^ etc.), for 
diflFerent VBO values for lens and disc type InAs/GaAs and InAs/InP quantum dots. 

For InAs/GaAs lens type quantum dot, the biexciton binding energy varies signifi- 
cantly and even changes its sign at VBO?^ 170 meV, leading to the unbound biexciton 
for VBO< 170 meV. Similarly for this system, binding energies of charged excitons vary 
over large range of values, with bound to unbound transition for X~ at VBO?^ 0.35 eV. 
For InAs/InP lens type quantum dot, in the "realistic"—"— range of InAs/InP VBO values 
(300 — 400 meV), charged complexes even reverse their relative position in a complicated 
way, e.g. for VBO= 0.3 eV, one observes the following ordering of spectral lines with the 
increasing energy: XX, X+, X~, X, while for VBO?^ 0.4 the order of lines is following: X~, 
XX, X+, X. For InAs/GaAs disc type quantum dot (" realistic"-^^--^^ YBO^ 200 - 300 meV), 
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one observes following order of spectral lines: X~, XX, X+ and X, while for InAs/InP disc 
type quantum dots (VBO^ 300 — 400 meV) one gets: X~, XX, X and X+. In all considered 
cases absolute binding energies and relative order of spectral lines depend significantly on 
the VBO, quantum dot shape and substrate material. And thus caution should be exercised 
when applying different VBO values in excitonic calculations. 

Binding energies must be calculated using the full many-body approach —^^i^: 
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where electron-electron (Jee), hole-hole {Jhh) and electron-hole J^h Coulomb integrals are 
calculated for electron and hole occupying the s-shell ground states, while the important 
correction due to correlation effects AEcorr can be attributed to the configuration mixing 
effects with higher lying states. In the Hartree-Fock approximation AEcorr = 0, while 
realistic values of binding energies {AEcorr 0) can be calculated using the full configuration 
interaction method^^^. 

For example, for lens type InAs/GaAs quantum dot, for the VBO= 0.1 eV, Jee = 
26.28 meV, Jeh = 19.54 meV and at the level of Hartree-Fock (perturbation theory) ap- 
proximation binding energy can be calculated as AE^^ {^~) = Jss ~ ^ss = 6.75 meV. 
This value is further reduced by the correlation correction AEf^^^ = 1.41 meV and there- 
fore finally E^~ = 5.34 meV. On the other hand, for VBO= 0.5 meV Jee = 20.88 meV, 
Jeh = 21.68 meV, leading to AE^^ {^~) = —0.8 meV. The CI calculated correlation cor- 
rection is AEf^^^ = 1.52 meV and thus we obtain E^ = —2.32 meV. 

The correction due to correlation mixing is important {^ 1.5 meV) and cannot be ne- 
glected, however its value does change significantly as VBO is varied. In the previous section 
I discussed that while the absolute values of excitonic fine structure splittings were results 
of the full configuration interaction procedure, the observed trends could be analyzed in 
terms of single exchange integral. Quite similar situation takes also place for the excitonic 
complexes binding energies. The absolute value of the binding energy must be calculated 
using the many body approach, yet its evolution with respect to the VBO can be under- 
stood at the level of Hartree-Fock (perturbation theory) approximation and several s-shell 
electron-hole Coulomb integrals (Jge, Jeh and Jhh)^ while the correlation correction AEcorr 
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is virtually unaffected by the choice of the VBO. 
V. SUMMARY 

Valence band offset is one of the important empirical parameters for the semi-empirical 
tight-binding method. By a thorough VBO analysis I have shown a non-trivial interplay 
between confinement potential due to the band offset and confinement due to strain. With 
strain effects artificially neglected, no shell structure for holes is present, with complicated 
charge density distribution due to light and heavy hole band mixing effects. With strain 
effects included, heavy holes and light holes are decoupled from each other by the biaxial 
strain and the shell-like structure characteristic for single band models is restored. Notable 
difference from this picture is observed for the strained InAs/GaAs lens type quantum dot. 
Due to anisotropy and inhomogeneity of strained induced confining potential splittings of 
excited hole states dominate over shell spacings, even for the largest considered VBO value. 
Peculiar structure of hole levels in InAs/GaAs in lens type QDs, so much different from 
InAs/InP QDs, origins thus from strain and not due to the different VBO as suggested 
by the other studies^^. On the other hand the characteristic shell structure of InAs/GaAs 
disc type QDs or InAs/InP lens type QDs also origins from strain effects. In all considered 
cases, even for "low-strain" InAs/InP systems, strain effects and large VBO value lead to 
spectral structure of hole levels known from experimental results. Therefore, strain cannot 
be neglected even in nominally "low-strain" InAs/InP quantum dots, whether disc or lens 
shaped. 

Valence band offset also affects many body properties such as excitonic fine structure and 
binding energies of excitonic complexes. Exciton fine structure splittings values reported in 
this work are order of magnitude larger than those reported by the EPM and thus much closer 
to the experimentally reported values without the necessity of inclusion of the "ordering" 
effects^. Fine structure splitting is also shown to vary depending on the matrix material 
(InP/GaAs) and the inclusion of strain. Charged excitons and biexciton binding energies 
vary strongly on the depth of confinement potential (VBO) with clear distinction from the 
simple effective mass ("harmonic oscillator") picture. This should be emphasized as the 
reliable prediction of excitonic binding energies is of key importance for so-called "inverse 
approaches"—. 
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The paper studies differences between flat (,, indium flushed") and lens type quantum 
dots, the relation between hydrostatic and biaxial strain in these systems, and quantitatively 
estimates the role of a jagged quantum dot/matrix material interface. The paper shows that 
the unrealistic p-shell symmetry of ground hole state would be obtained if ETB used the 
VBO value taken directly from the EPM and then discusses differences and similarities 
("strained band offset") between both approaches. Finally the paper shows that the VBO 
is not merely a technological parameter, and the caution should be exercised when applying 
different VBO values. 
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While the reported InAs/GaAs VBO extreme values vary from 60 meV to 500 meV, the 
"realistic" or the "recommended"^^ values can be limited to narrower range. 
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Lens QD (d=25 nm, h=3.5 nm). Strain effects included. 
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FIG. 1: (Color online) Single particle (electron and hole) energies (a) and effective gap Egap = 
ei — hi for InAs/GaAs lens type (d=25 nm, li=3.5 nm) quantum dot as a function of quantum 
dot (InAs) and matrix (GaAs) valence band offset (VBO). Single particle electron (b) and hole (c) 
energies are calculated with respect to electron (ei) and hole (hi) ground states energies. Patterned 
areas mark higher, excited states. 
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Disc type QD (d=16.8 nm, h=3 nm) 
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FIG. 2: Electron (conduction band - CB) and hole (valence band - VB) single particle energies 
calculated with respect to electron and hole ground states energies (ei and hi correspondingly) for 
InAs/InP and InAs/GaAs disc type (d=16.8 nm, h=3 nm) quantum dots as a function of quantum 
dot (InAs) and matrix (GaAs or InP) valence band offset (VBO). Strain-effects are either included 
or artificially neglected. 
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Lens type QD (d=25 nm, h=3.5 nm). 




FIG. 3: Electron (conduction band - CB) and hole (valence band - VB) single particle energies 
calculated with respect to electron and hole ground states energies (ei and hi correspondingly) for 
InAs/InP and InAs/GaAs lens type (d=25 nm, h=3.5 nm) quantum dots as a function of quantum 
dot (InAs) and matrix (GaAs or InP) valence band offset (VBO). Strain-effects are either included 
or artificially neglected. 
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Electron states of disc QD (d=16.8 nm, h=3 nm) 
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FIG. 4: Electron probability density isosurfaces for InAs/InP and InAs/GaAs disc type (d=16.8 
nm, li=3 nm) quantum dots as a function of quantum dot (InAs) and matrix (GaAs or InP) valence 
band offset (VBO). Strain-effects are either included or artificially neglected. 
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Disc type QD (d=16.8 nm, h=3 nm). Lens type QD (d=25 nm, h=3.5 nm). 
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FIG. 5: Splitting of the electron p-shell calculated for InAs quantum dots of different shape, 
matrix material and as a function of quantum dot (InAs) and matrix (GaAs or InP) valence band 
offset (VBO). 
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Electron states of lens type QD (d=25 nm, h=3.5 nm) 
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FIG. 6: Electron probability density isosurfaces in InAs/InP and InAs/GaAs lens type (d=25 nm, 
li=3.5 nm) quantum dots as a function of quantum dot (InAs) and matrix (GaAs or InP) valence 
band offset (VBO). Strain-effects are either included or artificially neglected. 
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Holes states of disc QD (d=16.8 nm, h=3 nm) 
Strain effects neglected 
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FIG. 7: Hole probability density isosurfaces in InAs/InP and InAs/GaAs disc type (d=16.8 nm, 
li=3 nm) quantum dots as a function of quantum dot (InAs) and matrix (GaAs or InP) valence 
band offset (VBO). Strain-effects are artificially neglected. 
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Hole states of lens type QD (d=25 nm, h=3.5 nm). Strain effects neglected. 
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FIG. 8: Hole probability density isosurfaces in InAs/InP and InAs/GaAs lens type (d=25 nm, 
li=3.5 nm) quantum dots as a function of quantum dot (InAs) and matrix (GaAs or InP) valence 
band offset (VBO). Strain-effects are artificially neglected. 
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Holes states of disc QD (d=16.8 nm, h=3 nm) 
Strain effects included 
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FIG. 9: Hole probability density isosurfaces in InAs/InP and InAs/GaAs disc type (d=16.8 nm, 
li=3 nm) quantum dots as a function of quantum dot (InAs) and matrix (GaAs or InP) valence 
band offset (VBO). Strain-effects are included. 



30 



Hole states of lens type QD (d=25 nm, h=3.5 nm). Strain effects included. 
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FIG. 10: Hole probability density isosurfaces in InAs/InP and InAs/GaAs lens type (d=25 nm, 
li=3.5 nm) quantum dots as a function of quantum dot (InAs) and matrix (GaAs or InP) valence 
band offset (VBO). Strain-effects are included. 



31 



> 

si 

>» 

o 
c 
liJ 



1.6- 



1.2- 



0.8- 



0.0- 



1 ■ 1 

n -□ -U ~U 

i 

/ 


( 






a =+1.0 eV 

V 


a =-1 .0 eV 

V 


VBO= 0.05 eV 


VBO= 0.21 eV" 


---- CB 


CB 


HH 


— HH 


----- LH 


— ^— LH 


O^O^O— 0— 0— o— o 




i^"~=A-=A— A— A— A— A 


V 







-5 5 

Position along growth axis (nm) 



FIG. 11: Strain-induced confining potentials for a InAs/GaAs lens type (d=25 nm, li=3.5 nm) 
quantum dot along [001] axis, calculated using the Bir-Pikus model and two sets of bulk VBO and 
ay parameters. Conduction band (CB) - squares, heavy-hole band (HH) - circles, light-hole band 
(LH) - triangles. 
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FIG. 12: Trace Tr (e) of tlie strain tensor and tlie biaxial B (e) = e^z — {^xx + ^yy) /2 component 
of strain for a disc type (d=16.8 nm, li=3 nm) and a lens type (d=25 nm, li=3.5 nm) InAs/GaAs 
quantum dot. Profiles has been calculated along [110] direction, z — 2 nm from the QD base. 
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FIG. 13: Strain-induced heavy-hole confining potentials for InAs/GaAs and InAs/InP disc type 
(d=16.8 nm, h=3 nm) and lens type (d=25 nm, h=3.5 nm) quantum dots. Profiles has been 
calculated using the Bir-Pikus model along [110] direction, z = 2 nm from the QD base. 
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FIG. 14: Excitonic fine structure for InAs/GaAs (open circles) and InAs/InP (squares) disc type 
(d=16.8 nm, h=3 nm) and lens type (d=25 nm, li=3.5 nm) quantum dots as a function of quantum 
dot (InAs) and matrix (GaAs or InP) valence band offset (VBO). 
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FIG. 15: Biexciton XX (black-circles), negatively charged X~ (red-squares) and positively charged 
X+ (blue-diamonds) excitons binding energies, calculated with respect to the single, neutral, ex- 
citon energy, for InAs/GaAs and InAs/InP disc type (d=16.8 nm, h=3 nm) and lens type (d=25 
nm, h=3.5 nm) quantum dots as a function of quantum dot (InAs) and matrix (GaAs or InP) 
valence band offset (VBO). 
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